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I. INTRODUCTION 

The gravitational interaction of n compact objects is 
a fundamental astrophysical problem. If one wants to 
tackle this problem within Einstein's general relativity 
[1], then one generally has to resort to numerical simu- 
lations, see, e.g., [2-4]. However, there exists a number 
of approximation methods. One of the most successful 
approximation schemes is the post-Newtonian (PN) ap- 
proximation, a slow motion and wide separation approx- 
imation. This allows an approximate solution of the field 
equations to some order, leaving only ordinary differen- 
tial equations for positions, momenta, and spins of the 
compact objects. It is convenient to encode these equa- 
tions of motion in terms of a Lagrangian potential or 
a Hamiltonian. We will calculate the PN approximate 
Hamiltonian via the canonical formalism of Arnowitt, 
Deser and Misner (ADM) [5]. 

In the present paper we concentrate on spin contribu- 
tions within the post-Newtonian approximation. More 
specifically we will derive the conservative n-body next- 
to-leading order (NLO) spin-orbit and spin(a)-spin(b) 
contributions, where a and b label different compact ob- 
jects, to the post-Newtonian Hamiltonian. These contri- 
butions were already derived for the binary case n = 2 in 
[6, 7]. Other derivations can be found in [8-12]. Next-to- 
lcading order spin-orbit contributions to the equations 
of motion for n — 2 were first obtained in [13] and 
essentially confirmed in [14]. The leading-order (LO) 
spin-orbit, spin(a)-spin(b), and spin(a)-spin(a) contribu- 
tions are well-known, see, e.g., [15-20]. For the next-to- 
lcading order binary spin(a)-spin(a) interaction see [21- 
25]. Some binary Hamiltonians of even higher order in 
spin can be found in [22, 26]. 

The spin contributions to the dynamics have to be sup- 
plemented by appropriate (i.e., sufficient within the ap- 



proximation scheme) point-mass 1 contributions. In gen- 
eral spin contributions can not directly be compared to 
point-mass contributions as the spin is a further expan- 
sion variable. However, for maximal spin, which is de- 
fined by a ratio of spin to mass-squared corresponding to 
the extreme Kerr solution, each power in spin is equiva- 
lent to half a post-Newtonian order. Let us recall that for 
maximal spins the leading order spin-orbit Hamiltonian is 
at 1.5PN order and the leading-order spin(a)-spin(b) one 
is at 2PN order, while formally counted (i.e., without re- 
lating the spin variables to the PN counting) both leading 
order Hamiltonians are at 1PN order. Similarly for max- 
imal spins the next-to-leading order Hamiltonians given 
in the present paper, formally having a post-Newtonian 
order of 2, are comparable to 2.5PN point-mass contri- 
butions for the spin-orbit case and to 3PN point-mass 
contributions for the spin(a)-spin(b) case. The point- 
mass dynamics to 2PN was completed for n = 3 in [28], 
for corrections see [29], and reduced to master integrals 
for arbitrary n in [30]. The spin-orbit Hamiltonian in 
the present paper thus completes the dynamics for max- 
imal spins and n = 3 to 2.5PN. (Notice that the dissi- 
pative 2.5PN point-mass dynamics can trivially be ex- 
tended from n = 2 to arbitrary many objects, see, e.g., 
[31].) However, it should be emphasized that the results 
in the present paper are valid for arbitrary n. Also notice 
that spins close to maximal are astrophysically realistic, 
see, e.g., [32]. 

Until the point-mass contributions to the post- 
Newtonian approximation arc not pushed to a higher 
number of objects n, the most useful application of the 
Hamiltonians given in the present paper is the investi- 
gation of the three-body problem with rapidly rotating 
objects in general relativity. This ideally fits to numeri- 
cal investigations as in [2-4], which are accurate beyond 
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1 In terms of covariant multipole moments [27], for a point-mass 
all multipoles except the monopole are neglected. The spin con- 
tributions in this paper arise from the covariant dipole moment. 
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the applicability of the post-Newtonian approximation 
but require much more computational power. An im- 
portant astrophysical application is the investigation of 
hierarchical triplets. The Hamiltonians provided in this 
paper allow an accurate treatment of, e.g., Kozai reso- 
nances [33, 34] in mergers of a black hole with a black 
hole binary when one or several of these black holes are 
rapidly rotating. One may also try to find stable solu- 
tions, such as the periodic ones for non-spinning objects 
given in [29, 35, 36]. Further the three-body problem is 
always interesting for the study of chaotic behavior. To 
foster such application the derived Hamiltonians for three 
compact objects are provided as Mathematica source files 
[37]. 

The paper is organized as follows. In Sec. II, we pro- 
vide a short introduction to the ADM formalism. After 
this, in Sec. Ill, a few details (namely constraint expan- 
sions, integration by parts and three-body integrals) of 
the calculation will be explained. In Sec. IV the results 
for the Hamiltonians and checks (and for the readers con- 
venience the appropriate center of mass vectors) will be 
provided. Last but not least, there will be some conclu- 
sions and further tasks given in Sec. V. 

The signature of spacetime is +2 in the present pa- 
per. Since the PN formalism is a perturbation theory 
around a flat Euclidean background it does not matter 
in principle whether the spatial indices in some tensor 
expressions are upper or lower ones (although the index 
position is important for the definition of some quanti- 
ties). These indices are denoted by small Latin letters 
from the middle of the alphabet (i, j,k, . . . ) and running 
from 1 to 3. Greek indices (/x, v, . . . ) are 4-dimensional 
indices running from to 3. Object labels are denoted 
by small Latin letters from the beginning of the alpha- 
bet (a, b, c, ...). In this paper we sum over all double 
indices (Einstein summation convention) except object 
indices. Sums over object labels are explicitly written in 
the expressions. Vectors are denoted by boldface letters 
and the scalar product between two vectors a and b is 
denoted by (ab) = (a • b). Furthermore, the speed of 
light c is set equal to 1 and if there is a (c _1 ) appear- 
ing it is just a bookkeeping parameter to get the correct 
post-Newtonian order of quantities or expressions. The 
Newton gravitational constant is denoted by G and we 
did not use a special convention for it. The reader may 
set G to a desired value in our expressions. 



II. THE ADM FORMALISM 

In the present paper, we will utilize the ADM canoni- 
cal formalism after gauge fixing [5], see also [38, 39]. At 
this stage the constraints of the gravitational field are 
solved (approximately in our case). The Hamiltonian is 
then given by the ADM energy expressed in terms of cer- 
tain canonical variables. The ADM formalism has shown 
to be valuable for calculating the conservative dynam- 
ics within the post-Newtonian and post-Minkowskian ap- 



proximations, sec, e.g., [40-42]. 

The constraints of the gravitational field can be written 

as 
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and arise as certain projections of the Einstein equations 
with respect to a timelike unit 4- vector with compo- 
nents ?v = (-N, 0,0,0) or n» = (1,-N l )/N. Here Ty 
is the induced three-dimensional metric of the hypersur- 
faces orthogonal to n^, 7 its determinant, R the three- 
dimensional Ricci scalar, Kij the extrinsic curvature, N 
the lapse function, N l the shift vector, y^jT^ v the stress- 
energy tensor density of the matter system, and ; denotes 
the three-dimensional covariant derivative. Partial coor- 
dinate derivatives di are also indicated by a comma. For 
the extrinsic curvature we used the ADM sign con- 
vention, i.e., 2NKij — — -fij.o + N i: j + Nj-j. 

In the ADM transverse-traceless (ADMTT) gauge de- 
fined by 



7T" = . 



(0) 
(7) 



which will be used throughout this paper, one has the 
decompositions 



in = ( 1 + i ) % + h V > 



n V = 7T^ T + jf« 



(8) 
(9) 



where hj^ and 7r^? T are symmetric and transverse- 
traceless, e.g, hj^ = hj^, hj^ = hj^j = 0. Notice 
that the form of the trace term in (8) is adapted to the 
Schwarzschild metric in isotropic coordinates, with ob- 
vious advantages for perturbative expansions. For our 
convenience we introduced a rescaled 4> = which is 
useful later in the expansion of the constraint equations. 
The longitudinal part, 7r y , of ir^ in Eq. (9) can be writ- 
ten in two equivalent forms, either in terms of tt 1 (which 
contains an inverse Laplacian A -1 ), 
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or in terms of V (which contains no inverse Laplacian) , 

(11) 
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The two vector potentials tt' 1 and V' 1 arc related by 



III. CALCULATION 
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and can be obtained as solutions of the momentum con- 
straint via 



(13) 



cf. Eq. (20). The transverse-traceless(TT) part of tt, tt^ t , 
is given by 



TtJ t — S M • i TT kl , 



with the partial space-coordinate derivatives di and the 
transverse-traceless projector 

= §[(*« - £sr 1 d i 8i){6 jh - &~ l djdk) 
+ (6a - A^didkXS^ - A" 1 ^) (15) 
-(Su-A^dMSij-A^didj)}. 

Now the four field constraints can be solved for the 
four variables and tt 1 in terms of hj^ , tt^ t and matter 
variables, which enter through the stress-energy tensor 
via the source terms % mattcr and %n iattcr . An analytic 
solution for cj> and tt 1 , however, can in general only be 
given in some approximation scheme. Finally, the ADM 
Hamiltonian -Hadm reads 
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This is the ADM energy expressed in terms of the canon- 
ical variables. The canonical matter variables are intro- 
duced in Sec. IIIB below. The canonical field variables 
are hj^ and tt^ t here, with the Poisson brackets 



{/^ T (x),4' T (x')} = 16nG5™5( X - x >) 



(17) 



Notice that beyond the post-Newtonian order considered 
here spin corrections to the canonical field momentum 
are needed [43, 44]. 



The field and source expansions starting at their lead- 
ing order are given by 
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where the subscript in round brackets denotes the 



(14) order. The hJ T field only occurs in leading order, namely 



(c ) . Since the TT field momentum is related to time 

derivatives of hj^ the leading order of 

The mass m a , canonical matter momentum P a , and spin 

variables S a are formally counted as m a ~ O (c _1 ) , 

P a ~ O (c -1 ) 3 , and S Q ~ O (c -1 ) 3 for dimensional rea- 
sons only (remember that for maximal spins one would 

have S a ~ O (c _1 ) instead). This counting comes from 
the fact that after setting c = G = 1 we require all quan- 
tities to be in units of length. Let us introduce symbols 
with a bar over them being the quantities in SI units and 
the other symbols the quantities in units of length, then 
it holds m a = ^fh a for the mass, t = ci for the time, 
Pa = f^Pa for the linear momentum, and similar for the 
spin variables. So the order counting comes from the c 
powers inserted to reconstruct the SI units. It should be 
noted that these counting rules will in general not give 
correct absolute orders in c if the SI units of the final 
expression are not taken into account. However, relative 
orders are always meaningful, which is all that is relevant 
for perturbativc expansions. Further notice that differ- 
ent counting rules are obtained if one assumes that all 
quantities are expressed in terms of mass units instead of 
length units when setting c = G = 1 , which is also often 
used in the literature. 



A. Constraint expansions 

The Hamilton constraint expansion is given by 
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where A = did{. To (c x ) 5 order, the momentum constraint can be expanded via 



(3) 'J 
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B. Source expansion 



The source of the field constraints % matter and % matter in terms of canonical variables were derived for spinning 
objects to linear order in the single spin variables and to the post-Newtonian order required in this paper in [45]. 
Higher post-Newtonian orders were treated in [44] and the formalism was worked out to all orders in [43]. The 
expansion of the source into powers of 1 /c reads 
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in the case of the Hamilton constraint sources and for the vectors as 
momentum constraint sources one obtains 
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Here P a i are the matter canonical momenta, £ a = (5(a; 1 — 
) with z l a the canonical position variable, and S a U)U\ = 
—S a (j)(i) is the canonical spin tensor. The latter is re- 
lated to the spin vector S a rp, by S a u\(j-\ = eijkS a (k) where 
Cijk is the Levi-Civita symbol. These variables have the 
canonical Poisson brackets 



sin v a cos <p a 
(S a (i)) = S a I sin 6 a sin <f> a 
cos 9„ 



(25) 



then possible canonical variables are the pairs <f> a and 
S a (3) = S a cos 6 a with 



{(f>a,S a(3) } = 1 



(26) 



all other zero, see [47] and also [6]. However, this intro- 
duces square roots as 
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(23) It is straightforward to check that (27) and (26) lead to 

(24) ^ 



all other zero. Notice that the spin- length 5 12 = S a (i)S a (i) 
is constant as all its Poisson brackets vanish. Therefore 
the spin has only two dynamical degrees of freedom. For 
some applications it is useful to work in a basis of phase 
space which makes this explicit, especially for investiga- 
tions regarding chaos [46]. If one parametrizes the spin 



C. Integration by parts 

The post-Newtonian expanded ADM Hamiltonian re- 
sults according to (16) from an integral over the right- 
hand side of (19). However, this integral can be greatly 
simplified. 
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First of all one can get rid of the ^'(3) 7r (5)TT term via in- 
tegration by parts, since is divergence free and one 



can rewrite tt 
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(5)TT 

in terms of derivatives of the V, 



(3) 



vector 



potential. Furthermore one can eliminate 0(6) "H^) via 
integration by parts and using Eq. (19c) in the system 
of constraint equations. One can eliminate the fr^ via 

rewriting ff^ in terms of Vfo derivatives as well and gets 

a source type term and another 0(2) (^(3)) contribution. 

After these integrations by parts one can change from 
a Hamiltonian in the TT degrees of freedom to a so called 
Routhian, a Hamiltonian in the particle degrees of free- 
dom and a Lagrangian in the propagating field degrees 
of freedom [40]. The TT degrees of freedom are then 
eliminated from the Routhian by inserting their approxi- 
mate solution. (The reason for this is that one can insert 
equations of motion for the time derivatives of the parti- 
cle variables appearing in the velocities of the TT degrees 
of freedom, which corresponds to a coordinate transfor- 
mation only [48].) The Hamiltonian resulting from this 
procedure is given by 
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where A -1 is the inverse Laplacian for usual boundary 
conditions. It is possible to rewrite Hjp N into a form 
where no point-mass part of hT^^ i s needed, see [45] for 
details. 



D. Three-body integrals 

Since there are at most three fields appearing (and 
no held which is generated by more than one body) in 
the integral for the Hamiltonian at 2PN spin-orbit and 
spin(a)-spin(b) level, namely 
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we will not get any integrals where a higher number of 
compact objects is involved. The abbreviations PP, SO, 
S, SS, and S 2 (some of them first appear later) stand for 
point-mass part (or point-particle part), spin-orbit part, 
spin part of the field, spin(a)-spin(b) part, and spin(a)- 
spin(a) part, respectively. (There is another three-body 
integral generated by two fields and one delta source in 
the integrand, ^^ a t yp C ' S °; for which we do not need the 
calculation procedures mentioned here and which is given 
in the results later.) So it is sufficient to calculate only 
three-body integrals for the n-body 2PN spin-orbit and 
spin(a)-spin(b) contribution to the Hamiltonian. We re- 
fer to integrals of the type mentioned above as three-body 
integrals, because they describe an interaction between 
three different position variables (if the object positions 
arc not distinct, one will get a two-body or one-body 
integral). 

The only three-body integrals appearing here are well- 
known and were already solved in three dimensions, 
namely 
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where s abc = r ab + r ac + r bc , r a = \x - x a \ and r ab = 
\i a —i b \. The first integral was solved in [49, Eq. (82,33)] 
and [50], and the second one in [40]. Note that the in- 
tegrals on the left-hand side of (36a) and (36b) are only 
formal expressions since they are divergent. Their (reg- 
ularized) solutions on the right-hand side are not unique 
and have to be fixed by certain consistency conditions, 
e.g., that Laplacians operating on different particle co- 
ordinates give certain functions. Further the integrals in 
the form given above are auxiliary functions, only their 
derivatives, which in fact are convergent, enter the phys- 
ical expressions. See discussion in, e.g., [40] for further 
details. 

After inserting these integrals, it was necessary to 
rewrite derivatives with respect to x into derivatives with 
respect to, e.g., i a to pull them out of the integral. 
(Derivatives with respect to components of particle co- 
ordinates z a are denoted by d\ a \) So in principle these 
(33) parts of the Hamiltonian can be calculated by integrating 



6 



appearing three-body integrals and afterwards differenti- 
ate them by the different particle coordinates three, four, 
or five times depending on the appropriate part of the 
Hamiltonian. The parts of the Hamiltonians which were 
calculated by the algorithm mentioned above are given 
by 
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Because of the two-body and one-body contributions in 
a three-body sum, one has to decompose the sum into 
a purely one-body part, a two-body part and a three- 
body part. The first two parts can be calculated by 
dimensional rcgularization procedures via the Riesz for- 
mula mentioned in [41] and the three-body integrals were 
calculated by applying the formulas given above. The 
delta-type one- and two-body contributions were calcu- 
lated by inserting Riesz kernel regulators and applying 
the method of dimensional regularization. The sources 
and the constraint decomposition in d dimensions will be 



given in a forthcoming publication. 

To eliminate the scalar products of n ab , n bc and n ac 
(here n ab = (z a — z b )/r ab ) one can make use of the fol- 
lowing identities 
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additional to n 2 ab = 1. Note the minus sign in front of the 
second expression, which comes from relabeling the first 
identity and changing the direction of the second factor 
in the scalar product. 

Furthermore one can remove one of the appearing unit 
vectors due to the fact that they are not linearly inde- 
pendent, namely 
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IV. RESULTS 

In this section we make use of xTensor [51], a free 
package for Mathematica [52] , especially of its fast in- 
dex canonicalizer based on the package xPerm [53], and 
some of our own code for evaluating integrals and deriva- 
tives. 

The Hamiltonians consist of several parts, which will 
be provided below, namely 





'ac 


rbc 


n a h = 


— n ac - 


— n bc 




Tab 


r a b 




Tab 


rbc 


n QC = 


— n ah + 


n 6c 










Tab 


, ^*ac 


n fcc = 


n ab 


+ — n 




rbc 





NLO 
SO 



[2],G 
NLO, SO 



H- 



[2],G 2 



rrNLO 
H SS 



-"nlo.ss 



NLO, SO 

[2],G 2 
NLO,SS 1 



[3],G 2 



H- 



NLO,SO 

[3],G 2 
NLO.SS ' 



(46a) 
(46b) 



where the numbers in square brackets denotes the num- 
ber of compact objects involved in the interaction. 



A. Spin-orbit Hamiltonian 

1. two-body interaction part 

The two-body interaction parts of the spin-orbit 
Hamiltonian linear in G and quadratic in G are given 
by 
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[2],G 
NLO,SO 



2],G 2 
NLO,SO 



E E *.2 



a b=£a 



' ab 



3_Pf_ _ 3 (n Qb P b ) 2 _ 5mfcP| _ 3 (P a P fc ) _ 3 (n ab P a )(n ab P b ) 

4 m a m b 2 m a m b 8 4 4 m 2 

((n ab x P b )S a ) 



((n ab x P a )S a ) 



EE 

a b=£a 



(P a P b ) | 3 (n ab P a )(n a& P b ) 

m a m b m a m b 
G 2 



3 (n ab P a ) _ 2 (n afc P fc ) 

4 to?, m a m b 



((P„ x P b )S a ) 



ab 



11 TO \ / 15 

y?7i fc + 5^- j ((n ab x P a ) S a ) + ( 6m a + y?7i b ) ((n ab x P b ) S Q ) 



, (47) 
(48) 



and consist of two-body interaction parts of the field integral H^)^^) 2 an d a part of the delta-type integrals. 
Comparison with l/r 2 2 and 1/Vi 2 terms of the next-to- leading order spin-orbit Hamiltonian in [6, 45] leads to full 
agreement. 



2. three-body interaction part 
The three-body interaction Hamiltonian is always at G 2 level and consists of three different parts: 



H 



[3],G 2 



rr mattcr,SO , rr mattcr.SO , rrTT-part,SO 



'NLO,SO — JjI o-typc " 0(2)tt(3) 2 ^ - , - , a</>(2) 2 /i TT (4) • 

The delta-type part (which results from the source parts of the integrand) is given by 



(49) 



H 



matter, SO 
<5-type 



G ( 5 1 \ m b m c 



EE E t- — + - 



a b^a c^a,b 



((n o6 x P a )S a ) 



(50) 



The pure field part is given by 



H 



matter.SO \ " \ " \ ' r <2^ 
0(2)tt(3) 2 - U TO ° 

a &7^a c^a,b 



G 



8 8 



4r ab 4r Q 



2 2 2 

^abc ^ \ ^*^ c Tab^ac Tab^bc ^ac^bc ^*acT^ c ^ab^bc 



((n bc x P b ) S c ) 



^"ab^ac ^"ab^bc ^"ac^bc ' ' ac^bc ^ab^ ac 



3^ab 3rfc \ - ) ((n ac x n bc ) S c ) f 16(n ae P b ) 

- — + — t ) (( n - x p w s =) r + 3 ( — — — 



12 2_ ^ _5 2r ab ^ 7rv _ 2r ab _6r« 

.2 , 



1 ab 



, r ao r ac r bc r QC ■ ac 

The TT field part of this Hamiltonian is given by 



r bc 




r 2 
'ab 


T ab r ac 




Sr h ) 


- + 

oc 


r *br 2 a J 



abc 

(n ab P b ) 
(n bc P b 



(51) 



H 



TT-part,SO _ 
3</>(2) 2 /i TT (4) — 







--('< 

Tab \ 


1 r bc \ 




l b abc y 


T ac 


' 2r a J 


1 ab X 



q% m a m b 

2-^2-^2^ m c 

a b=£a c=£a,b 

+ 4" (- 2r lc + r 2 bc - \r ac r bc + ^ ) }>(n ac P c )((n QC x n bc ) S c ) 

' ~u \ £ £ Tac , 



' be 



' ab 

J_J_L 



i 



Tab 



Sabc \ 8rl c 4r ac r bc ^r 2 ac r bc r ab \ 8r ac r bc 8 r, 



5 3 r ac 1 r bc 



, ab v 8 4r bc 8r 2 ac 
+(a o 6) 



1/3 1 , 3r bl 

<7.i 

1 /l 



06 



3 f ac _ \ 1 \c_ _ 1 r bc 

4 r bc 8 r QC 8 r 2 ac 



((n ac x P c ) S c ) 



(52) 



where (a <-> b) denotes an exchange of object labels a and b in all of the preceding terms. 
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B. Spin(a)-Spin(b) Hamiltonian 

1. two-body interaction part 

The whole spin(a)-spin(b) Hamiltonian linear in G (which is a sum of the delta-type part, the TT-part, and the 
vector potential part) is given by 



H- 



2],G 
NLO,SS 



1 



Am a m b 



6((n ab x P b ) S a )((n ab x P a ) S 6 ) + -((n ab x P a ) S a )((n ab x P b ) S b ) 



a b^a ab 

-15(n Qb P Q )(n ab P b )(n ab S Q )(n ab S b ) - 3(P a P b )(n Qb S a )(n ab S b ) + 3(n ab P b )(P Q S Q )(n ab S b ) 
+3(n ab P a )(P b S a )(n ab S b ) + 3(n ab P b )(n ab S a )(P a S b ) + 3(n ab P a )(n ab S a )(P b S b ) 

-i(P b S a )(P a S b ) + (P Q S a )(P b S b ) - 3(n Qb P Q )(n Qb P b )(S a S b ) + i(P a P b )(S a S b ) 

-((n o6 x P a ) S a )((n ab x P a ) S 6 ) - (n ab P a )(n ab S a )(P Q S b ) + (n ab P a ) 2 (S a S b ) 



2m 2 



The G 2 two-body interaction Hamiltonian is given by 



H 



[2],G 2 
NLO,SS 



EE 

a b^a 



6G 2 m a r 



ab 



(S a S b )-2(n Qb S a )(n ab S b ) 



(53) 



(54) 



We neglected appearing S 2 terms due to consistency reasons (the stress-energy tensor does not contain S 2 expressions) 
Note that these two parts of the Hamiltonian are also in perfect agreement with [7, 12, 45]. 



2. three-body interaction part 



The three-body interaction part of this Hamiltonian is given by ss = H^^*f£f, namely 



H 



matter, SS 
*(2M3) 2 



EEE 

a b^a c^a,b 

2 



abc 



r 



ab 

i 

' ab 



— 2r ac + SS. 

Tbc 



((n x n bc ) h 

a, )((n ac 

x n bc ) S 6 ) 
2 



'ac'bc ^ab^ac 



4 

r 2 " 

' ab 



(n ac S a )(n 

ac^b 



'ac'bc 



2r„ 



2r bc 



+ (n bc S a )(n ac S b ) 



r 



ac 

He 
1 

~2~ ' 



6 



1 ab 



2rt 



r 



1 

Tab 

'ac'bc 
2 









\ ' ac 




4( 



1 - 



2r, 



r bc 

ac 



'V 



TV 



ac 

r bc 

1 tic 



~i- 1 



1 ao 



2r 



ac 
3 

a c 
2 
be 



'ac'bc 

, 3(r a 



2r ab 

+ r bc ) 2 3(r ac + r bc ) a 



2r2 



ac' oc 



3 

r ab 

2 

r 2 " 

' ab 



^ r ac + 3r ac r bc + r ac r bc 



ac 
Tbc 



' ac 



3(r 



) 2 , 3(r ac + r bc ) 3 



2^ b 



' a6 



V 



ab 



r ac r bc 



9 

-r, 



4r 3 r? 

^' ac' 6c 
,2 



_ 

?"bc 



Z 'ab 



?~ab 



(S Q S b 









\ ' ac 




'ab \ 



3r ac 



3 r ac 

2 r 3 " 
z 'be 



r bc 



■ f ab 
3 

ac: 

.3 

be 



be 

2 

r 2 " 

1 



+ (n QC S a )(nb c Sb)<j + 



2r 3 



ab 



r ac r bc 



ac ac ac \ I , / . . i\ 



2r 3 

z 'be 



(55) 



C. Approximate Poincare algebra 

As a check of the Hamiltonians given in the previous section we look at the global Poincare algebra, sec, e.g., [6]. 
For this we need the center of mass vector 

G = -^/d»xxA^. (56) 
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Since the center of mass vector integrals arc given by the Hamilton constraint equations which are one order below 
the appropriate integrals for the Hamiltonians, there are no explicit three-body parts appearing there. For the readers 
convenience we provide them here as well. (Notice the abuse of vocabulary, in fact G/H is the center of mass, but we 
refer to G as center of mass vector.) The Newtonian center of mass vector can be calculated trivially and is given by 



E m a z a 



(57) 



The 1PN point-mass center of mass vector is given by 

>2 

a 

2m a 



G 



1PN 
PP 



- E 9™" Za E E 



Gm a m b . 



2 '—' f— ' r ab 

a b=£a 



see, e.g., [54, 55]. The leading order spin-orbit center of mass vector is given by 



G s8 -E^r( p - xS ») 



(58) 



(59) 



There exists no leading order spin(a)-spin(b) center of mass vector. The next-to-leading order spin-orbit center of 
mass vector is given by 



/~iNLO 



a b=ta 



r a b ^m a 



-5(P a xS Q ) + ((P a xS Q )n ab ) 



5z a + z b 



Tab 



EE 



G_ 

Tab 



^(P fc x S ) - \{-a.ab x S Q )(P b n ab ) - ((P 6 x S a ) n ab ) Z ° + Zfc 
2 2 r ab 



a b^ta 

and the next-to-leading order spin(a)-spin(b) center of mass vector is given by 



^ss 



G 



EE 

a b^a 



S Z 

(S b n ab )-^- + (3(S a n ab )(S b n ab ) - (S S&))-#- 



ab 



ab 



(60) 



(61) 



see, e.g., [6, 7, 45] (notice that there is a misprint in Ggo° in the published version of [45]). For comparison of the 
G° parts (up to linear order in spin) of the center of mass vectors and Hamiltonians one can use the center of mass 
vector and the Hamiltonian which can be calculated by integrating the source % matter directly and setting 7^ = Sij, 
which results in 



G 



SRT 



E 



P a eS a (i)(£) 



i a + jml + PI 



ffSRT = E VV, + • 



(62) 
(63) 



From the references given above and, e.g., [29, 54, 56] one can also get the Hamiltonians needed, which are given by 



= 



1PN 
PP 



ttLO 



E 



1 y-v y-, Gm a m b 



E 



2m a 2 j—* r ab 



G_ 

Tab 



(Pi) 2 

8m3 



EE E 



EE 

a b^a 

G 2 m a mbm, 



f^ps + i (7(P a P b ) + (n ab P a )(n ab P b )) 
2m a 4 



ttLO _ 



a b^a c^a,b 

3m b 
2m, 



E E r 2. 



^^ab^ac 

((n ab xP a )S a )-2((n ab xP b )S a ) 



a b^a '<**> ^ 

_C 

2r 



E E 2r 3 

a b^a ab 



3(n ab S a )(n ab S b ) — (S a S b 



EE 

a b=£a 



G 2 m 2 a mb 
2r 2 



(64) 



(65) 
(66) 
(67) 



r 



It is now straightforward, though rather lengthy, to check for checking spin-orbit parts of the {G, H} part of the 
that the global Poincare algebra is fulfilled. Note that 
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Poincare algebra for example one also has to include 
{Gso,-ffso} due to the spin Poisson brackets. 

The Hamiltonians given above are necessary to check 
the Poincare algebra relations involving the derived 
Hamiltonians. They are not sufficient to simulate the 
full post-Newtonian dynamics at 2.5PN. Additionally to 
the derived next-to-leading order spin-orbit Hamiltonian 
and the Hamiltonians mentioned above, one needs the 
three-body 2PN point-mass Hamiltonian [28, Eq. (5)] 
and [29, Eq. (A.l)], the leading order spin(a)-spin(a) 
Hamiltonian 



ttLO 

H S 2 



G mb 



a b^a ab 



C c 



3(S a n afc ) 



(68) 



given in, e.g., [25, Eq. (13)] (the constant Cq q 
parametrizing the quadrupole deformation due to spin 
for the ath object, with Cq a = 1 for a black hole), and 
the radiative 2.5PN point-mass Hamiltonian provided in 
[31, Eq. (41)]. 



D. Another derivation of the Hamiltonians 

Because of the momentum independence of the three- 
body part of the spin(a)-spin(b) Hamiltonian, one cannot 
check this part using the Poincare algebra. Therefore 
we rederived the Hamiltonians H^q = Yl a ^a(4)S a and 
_ff|g LO using the formalism given in [6] via the precession 
frequency ^(4)' ^1- (4-10) m [6]; and compared this with 
our result given above. (Note that for the spin(a)-spin(b) 
Hamiltonian it is in principle not necessary to derive all 
f2 a . One only needs the term containing one of the spins. 
So one has to multiply a factor 1/2 when adding up all 
parts of the Hamiltonian namely = \ J2 a ^(4)aS a 

to avoid overcounting.) Both results for spin(a)-spin(b) 
and spin-orbit are identical with our previous results. It 
is explained in the appendix why they should not even 
differ by a canonical transformation. The field variables 
necessary for calculating the spin precession frequency 
were taken from [40, 45, 57]. 2 

Notice that the performed rederivation of the Hamil- 
tonians via the spin precession frequency provides a very 
strong check because it needs expressions for lapse and 
shift, which are eliminated in the formalism we used be- 
fore. 



V. CONCLUSIONS AND OUTLOOK 

We have derived the post-Newtonian next-to-leading 
order conservative spin-orbit and spin(a)-spin(b) gravita- 



tional interaction Hamiltonians for arbitrary many com- 
pact objects. The spin-orbit Hamiltonian completes the 
knowledge of Hamiltonians up to and including 2.5PN for 
three compact rapidly rotating objects. The Hamiltoni- 
ans were checked with the help of the Poincare algebra 
and rederived with the independent method from [6]. 

A possible astro-physical application of our computa- 
tion should be the exploration of Kozai resonance in hi- 
erarchical triples containing spinning compact objects in 
a fully post-Newtonian accurate manner. Recall that in 
hierarchical triples experiencing Kozai resonance the or- 
bital eccentricity of the inner binary secularly evolves, 
mainly due to the tidal torquing between the inclined 
inner and outer orbits [33]. And, the general relativistic 
periastron advance of the inner binary can interfere with, 
and in principle terminate, the evolution of its eccen- 
tricity [34]. Therefore, the present computation should 
be useful in extending the detailed analysis presented in 
[34]. Interestingly, we note that Kozai resonance, as dis- 
cussed in Ref. [34] , is also proposed as a scenario to merge 
massive binary black holes resulting from galaxy mergers 
[58]. 

To best of our knowledge not even the leading or- 
der spin Hamiltonians were applied in this context. 
Though the next-to-leading order effects derived here 
are considerably weaker within the validity of the post- 
Newtonian approximation, they can still be important. 
For the three-body case many configurations are poten- 
tially chaotic and weak interaction terms can have big 
effects. Further the leading order spin Hamiltonians only 
consist of two-body interactions, i.e., the objects interact 
pairwise with each other as in Newtonian gravity. At the 
next-to-leading order the most complicated parts of the 
Hamiltonians are three-body interactions and thus pro- 
vide not just a refinement of the leading order dynamics. 
At the next-to-lcading order the complexity of Einstein's 
theory of gravitation becomes apparent. Finally, the size 
of next-to-leading order effects provides a handle on the 
accuracy of the leading order. 

Additionally, it is interesting to see whether one can ex- 
tend known three-body solutions without spin, see, e.g., 
[29, 35, 36], to the three-body problem with spin at cer- 
tain order. In the literature, there exist paramctrizations 
for the binary case at leading order spin-orbit [59, 60]. 
For special configurations like spins aligned to orbital an- 
gular momentum, a parametrization for three bodies in- 
cluding next-to-leading order spin-orbit interaction seems 
to be possible, see [61] for the binary case. 

To foster application of the derived Hamiltonians we 
provide them for three compact objects as Mathematica 
source files [37]. 



Note there is a misprint in Eq. (5.6b) in [57] in the term 
2 — - — {. . . }. This term should be 2—^— {. . . } and in the last 



term before I^p the free index of nj, should be i instead of j. 



Appendix A: Relation to the spin variable used by 
Damour, Jaranowski, and Schafer 

In the past [45] the method in [6] was already used as 
a check, but it was not clear why the Hamiltonians were 
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in perfect agreement (and not differing by a canonical 
transformation). To explain this issue we compare the 
canonical spin S^ JS used in [6] and the canonical spin S a 
used in the present paper and in [45]. 

The comparison was done in the following way. We 
constructed the matrix due to Eq. (2.7) in [6]. From 
that we calculated the symmetric matrix square root H lJ 
which relates S„ JS to the spatial components of the co- 
variant spin 4- vector S a u (which fulfills the covariant spin 
supplementary condition S ali u% = 0). 3 Notice that it is 
enough to compare the definitions of the canonical spin 
variables, as the formalism in [6] is based on the spin 
equation of motion, in which corrections to the canonical 
position and momentum are of higher order in spin and 
can be neglected. 

Now we split up the equation which relates S a u to the 
spin tensor S£ u given by 



Sail — 2 v 9(A) e pvaf}u' / a S" , 

in a (3+1) manner. This gives 



(Al) 



&ai = UjkT 7 



Pa m (yiSaJr 



2 flPa Sa i 



(A2) 



using the (3+1) decomposition of n M , e 0123 = 1 such that 
£oi23 = -1 and so e oij k = -eyjfc, and u v a = P^lm a . Af- 
ter the (3 +1) split we insert (nS a ) l = -P a k"f kj S ajl /m a , 
nP a = — \J m\ + 7^ P a iP a j and the transformation from 
the covariant spin to the Newton- Wigner spin, namely 

c _ C Pai{ n S a )j Pgj(nS a )i 

O a ij — Daij — —f^- -f — —f^- ■ (AO J 



nP a 



nP n 



Now one has to go from the Newton- Wigner spin tensor 
in a coordinate basis S a ij to the canonical spin tensor in 



a triad basis S, 



via Sat 



. This 



canonical spin tensor can be related to the spin vector 



S a via S a 



or S, 



a(i) — 2 ei i k Sa(])(k) 



The transformation going from S a to S a using (A2), 
(A3), the basis transformation, and the relation be- 
tween spin tensor and spin vector can be compared with 
(H^ 1 )' 1 ! calculated perturbatively from H 1 * mentioned 
above. From this calculation one can see that there is a 
deviation from the canonical spin used in [6] and S used 
here of the form 



=.DJS 
ai 



£»«) + 



TT 



Pa j (PaS a ) 



P /i TT P . . 

Jr aii<'(4) mn - r a m^a (n) 



(A4) 



where the appearing field variables on the right-hand 
side have to be evaluated and regularized at the posi- 
tion z a of the ath object. This deviation is three post- 
Newtonian orders after the leading order. Note that the 
lengths of the spins have to be equal (which one can 
see from (A4) is fulfilled at the order considered), since 
(S° JS ) 2 = Sl = S 2 a with 2S 2 a = S^S auv . The deviation 
in (A4) is beyond the post-Newtonian order considered 
in the present paper, thus showing why the Hamiltoni- 
ans calculated in the ADM formalism are exactly iden- 
tical to the Hamiltonians calculated via spin precession 
frequency. 
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